NT@UW-09-22 



Isospin-symmetry-breaking corrections to superallowed Fermi P decay: 

Radial excitations 
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Based on an exact formalism, we study the effects of isospin-symmetry-breaking interactions on 
superallowed 0^ 0"^ transitions. We calculate the second-order renormalization of the Fermi ma- 
trix element due to radial contributions and show that radial excitations neglected in the treatment 
of Towner and Hardy are significant. These are estimated to decrease the isospin-symmetry-breaking 
corrections. Our results provide a correction term that can be included in existing approaches. 
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I. INTRODUCTION 

Superallowed nuclear (3 decays play a key role for pre- 
cision tests of fundamental symmetries: They provide 
the most stringent test of the conserved-vector-current 
(CVC) hypothesis, the best limit on scalar interactions, 
and the most precise value for the up-down Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Vud [H- The 
precise extraction of Vud from superallowed transitions 
requires the evaluation of theoretical corrections due to 
isospin-symmetry breaking (ISB) 0, Q and radiative Q 
effects with an uncertainty of 10%, to guarantee a de- 
sired accuracy of 0.1% for Vud- Because of the very high 
precision reached experimentally, the uncertainty of Vud 
is currently limited by ISB and radiative corrections [l|. 
This presents a challenge for nuclear theory. 

Superallowed 0"*" — > 0"*" Fermi (3 decays depend only on 
the vector part of weak interactions, and with CVC, the 
transition "/t value" should be nucleus-independent: 



ft 



2n^h'^ In 2 



const. 



(1) 



where Gy is the vector coupling constant and Mp is the 
Fermi matrix element. The CVC hypothesis depends on 
the assumption of isospin symmetry, which is not exact 
in nuclei, but broken by electromagnetic and quark mass 
effects. As a result, Alp is reduced from the symmetry 
value, Mq = V2 for T = 1 parent and daughter states, 
by ISB corrections Sc, 



\Mf\'^\Mo\'{1-Sc). 



(2) 



In addition, there are radiative corrections, but we focus 
on Sc in this paper. 

Towner and Hardy (TH) have shown that the cal- 
culated ISB corrections eliminate much of the consider- 
able scatter present in the uncorrected ft values, and the 
statistical consistency among the corrected J^t values is 
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some evidence that the calculated ISB corrections are not 
unreasonable. However, the precise extraction of Vad and 
the importance of testing the Standard Model have stim- 
ulated us to undertake a reevaluation based on an exact 
formalism Q . In addition, our goal is to connect the cal- 
culated corrections to the accurate understanding of ISB 
in nuclear forces [1, 0| • Reference 0] showed that there 
are specific corrections to the treatment of TH and 
explained these using schematic models without making 
numerical estimates. In particular, radial excitations are 
neglected by TH, and as a result the transition opera- 
tor violates the isospin commuation relations. The TH 
approach is discussed in more detail in Sect. [Ill 

Following our work, there have been a number of 
theoretical developments: Auerbach studied ISB correc- 
tions assuming that the dominant physics is due to the 
isovector monopole resonance @. Assuming that cer- 
tain reduced matrix elements of the isosvector part of 
the Coulomb potential are identical, Auerbach showed 
that ISB corrections vanish unless one takes into account 
the energy differences between components of the isovec- 
tor monopole state of different isospin. Using schematic 
models, the resulting estimates for Sc are considerably 
smaller than the TH results. Liang et al. carried out 
random-phase-approximation calculations based on rel- 
ativistic density functionals. They also obtain smaller 
values for Sc and their corrected Tt values are statisti- 
cally consistent In addition, Satula et al. have ana- 
lyzed isospin mixing and implemented isospin projection 
in density-functional calculations [loj . 

At the same time, the experimental precision has been 
improved in several cas es llj and a recent branching ratio 
measurement for '^^Ar [Il[ has improved the ft value for 
this T = 2, = — 2 decay to 0.8%, a precision that is 
ncaring the 0.3% standard for inclusion into the set of 
well-known T = 1 decays. This is especially interesting, 
because ISB corrections appear to be larger for T = 2 
superallowed transitions. 

In this work, we start from the exact formalism devel- 
oped in Ref. which is reviewed in Sect. lIIII In Sect. II VI 
we write the correct isospin operator as a sum of the TH 
operator and a correction term due to radial excitations. 
We show that the neglected radial contributions are sig- 
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nificant at leading (second) order in ISB interactions. In 
addition, for certain conditions the correction term can- 
cels ISB corrections to Fermi matrix elements. This is 
similar to the case of Auerbach when the components 
of the isovector monopole state are degenerate [8[ . Our 
results demonstrate explicitly the importance of radial 
excitations and wc discuss the implications in Sect. fVl 

II. TH APPROACH TO ISB CORRECTIONS 

With CVC, the matrix elements of weak vector interac- 
tions in nuclei are not modified by nuclear forces, except 
for corrections due to ISB effects. Therefore, one has 
to evaluate the contributions from electromagnetic and 
charge-dependent strong interactions to the Fermi ma- 
trix element Mp = {f\'r+\i) between the initial and final 
states for superallowed (3 decay, \i) and |/), respectively. 
Here t+ denotes the isospin raising operator. 

Towner and Hardy use a second-quantized formu- 
lation to write the Fermi matrix element as 

MF=Y^{f\aiap\i){a\T+\P), (3) 

where a]^ creates a neutron in state a and ap annihilates 
a proton in state (3. Thus, the label a is used to denote 
neutron creation and annihilation operators, while (3 is 
used for those of the proton. This notation is different 
from the standard notation in which ha is used to 
denote proton annihilation operators. The single-particle 
matrix element (a|r_|_|/3) is assumed to be given by 

/•oo 

{a\T+\(3) ^ 5a,0 Rl{r)Rl{r)r^dr = 5a,pra, (4) 
Jo 

where i?"(r) and R'^ir) are the neutron and proton radial 
wave functions, respectively. Because the radial quantum 
numbers of the states a and /3 are set to be the same, 
Eq. (ID) assumes that r+ creates a neutron in the single- 
particle state with the same quantum numbers as those 
of the annihilated proton. This is not the case in the 
presence of ISB. As a result, the operator in Eq. ^ is 
not the correct isospin operator and the Standard Model 
isospin commutation relations are lost. 

We observe that Eqs. ([3]) and ([4]) correspond to the 
second-quantized isospin operators 

a, 13 

T- = tI^^ Sa,i3 rl a!^paa , (6) 

a, 13 

SO that the commutation relations are given by 

[t+,t^] =^|rapaJ,aa-^|r/3pa]3a/3 7^To, (7) 

a (3 

which shows explicitly that the Standard Model isospin 
commutation relations are violated if one uses the isospin 



operator of TH. The violation is due to neglecting parts 
of the proton wave functions that are in radial excitations 
when expanded in a neutron basis. This formal problem 
motivated us to develop an exact formalism for ISB cor- 
rections and to study corrections to the TH treatment ■ 
For the calculation of 5c, TH proceed by introduc- 
ing into Eq. ([3]) a complete set of states for the {A — 1)- 
particle system, |7r), which leads to 

MF^Y.^f\aiU){Aac.V)rl. (8) 

a,7r 

The TH model thus allows for a dependence of the radial 
integrals on the intermediate state tt. If isospin were an 
exact symmetry, the matrix elements of the creation and 
annihilation operators would be related by hermiticity, 
{T:\aa\i) = (/laalTr)*, and all radial integrals would be 
unity. Hence the symmetry-limit matrix element in the 
TH model is given by 

M, = Y,\{fwi\^)?- (9) 

a.TT 

Towner and Hardy divide the contributions from ISB 
into two terms. First, the hermiticity of the matrix ele- 
ments of Cq, and a]^ will be broken, and second, the radial 
integrals will differ from unity. Assuming both effects are 
small, TH separate the resulting ISB corrections into two 
model-dependent parts 01 

Sc = 5ci + 5c2 , (10) 

where in evaluating 5ci all radial integrals are set to unity 
but the matrix elements are not assumed to be related 
by hermiticity, and in evaluating 5c2 it is assumed that 
(^|ac|*) = but rl ^ 1. 

While Eq. ([3]) accounts for important effects of the 
Coulomb interaction on the radial wave functions, ra- 
dial excitations are neglected in Eq. We will show 
in Sect. llVl that radial excitations contribute at the same 
(second) order in ISB interactions and are estimated to 
decrease the ISB corrections 5c2 obtained using the TH 
isospin operator. 

III. EXACT FORMALISM AND THEOREMS 
FOR ISB CORRECTIONS 

In this section, we recall aspects of the formalism and 
theorems of Ref. @ . These show that there are no first- 
order ISB corrections to Fermi matrix elements and pro- 
vide a basis for a complete second-order calculation. The 
formalism starts from the correct isospin operator 

r+=^al6„, (11) 

a 

where a represents any single-particle basis, and aj^, cre- 
ates neutrons and ba annihilates protons in state a. The 
Fermi matrix element is then given by 

MF = {f\T+\i), (12) 
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where |i) and |/) are the exact initial and final eigenstates 
of the full Hamiltonian H = Hq + Vq, with energy Ei 
and Ef, respectively. Here Vc denotes the sum of all 
interactions that do not commute with the vector isospin 

= [Vc, T] ^ and [Hq, T] = . (13) 



operator T : 
[H,T] 



We use round states, [i) and |/), to denote the bare ini- 
tial and final eigenstates of Hq. Obtaining these states 
requires a solution of the A-body problem. 

The full initial and final eigenstate |i) and |/) can then 
be written as 



I/) 



I/) 



1 



E, - A^HA, 



1 



AfHA 



(14) 



(15) 



with projectors = 1 - [i){i[ and A/ = 1 - |/)(/|. The 
unperturbed states |i) and |/) are related by an isospin 
rotation, so that |/) is the isobaric analog state of |i): 

|/) = -Lr+|z) and |z) = -Lr_|/), (16) 



2T 



2T 



where T denotes the isospin of the states |i) and |/), 
and for simplicity we consider parent cases with isospin 
projection T. = (iV - Z)/2 = -T or T = 1, 0, so 

that Ma = (/|t+|«) = VW. The factors Zi and Zf are 
taken to be real and ensure that the full eigenstates are 
normalized. With (/| t+A^ = and Ajt+ |«) = 0, the 
exact Fermi matrix element is given by 



A/o + (/|FcA/ 



1 



Ef -AfHA 



1 



E, - A,i/A. 



(17) 



This is the first theorem Q and, because of Z^./ = 1 + 
0{Vq), demonstrates that there are no first-order ISB 
corrections to Mp. 

Another formulation is obtained by expanding in the 
difference of the charge-dependent interactions AVc be- 
tween the initial proton-rich and final neutron-rich states. 
In this case, the full Hamiltonian is given by 



H = Ho + AVc, 



(18) 



where Hq includes the effects of Vc common to the ini- 
tial and final states, for example the Coulomb interac- 
tions in the core, and AVc takes into account all charge- 
dependent interactions of the extra proton with the other 
nucleons in the initial state. In this formulation, the bare 
states, [i) and |/), are eigenstates of Hq, but are not 
eigenstates of the isospin operator: 

K) = E 7T'|T',T,), (19) 

T'^\T^\ 

I/) = E 7T'|T',T, + 1). (20) 

T'^\T^ + l\ 



In the isospin-symmctry limit, jt' = ^T'.t, where T is 
the isospin of the bare states of Eq. (which is also 
the dominant isospin in the presence of ISB). 

In this case, the full eigenstates can be written as 



|/) = |/) and [i)^^Z[i) + 



1 



E, - A.HqA, 



■^,AVc\i) 



(21) 



and we obtain for the exact Fermi matrix element 



Mf = E l7T'lVr'(T' + l)-T,(r, + l) 

T' 

1 



(/|r+A, 



Ei — AiHoAi 



A,AVc[i), (22) 



where the sum is over T' ^ max{[T^[,[T^ + 1[). The 
|T, Tz) expansion of the states [i) and |/) presents a more 
careful evaluation of the second theorem of Ref. Q . Since 
7T = 1 + 0{V§) and -fT'^T = 0{Vc), it follows that 
(/|r-|_Ai is of first order in ISB interactions. Combined 
with Z = 1 + 0{{AVcf), Eq. ^ explicitly shows that 
ISB corrections to Mp start at second order. In the fol- 
lowing, we will work with the first formulation. 



IV. RELATION BETWEEN TH OPERATOR 
AND ISOSPIN 

Next we derive a relation between the TH operator and 
the correct isospin operator based on the exact formal- 
ism. We use basis states given by conveniently-chosen 
one-body potentials U and U + Uc , where Uc accounts 
for charge-dependent effects. The single-particle (sp) po- 
tentials are chosen to minimize the effects of residual in- 
teractions, 



AV = V + Vc -{U + Uc), 
so that the Hamiltonian is given by 



(23) 



H = Hsp + AV with Hsp^T + U + Uc- (24) 

We express the isospin raising operator in a mixed 
representation, where [a) denotes the eigenstates of the 
single-particle Hamiltonian Hgp and [a) the eigenstates of 
the isospin-symmetric part T+U. The creation operators 
in the two bases are related by 



E 4' ("'I") 



(25) 



where the tilde indicates the basis and the sum is over 
all single-particle quantum numbers. The correct isospin 
operator, Eq. pT|) . can then be expressed as 



= E] "^^z {a' [a) ba ■ 



(26) 



a, a' 



which TH use as their starting point. However, Eq. 
allows the states |a) and |a') to have different radial 
quantum numbers n and n', to be explicit 



nljm and a' — n'ljm , 



(27) 
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with orbital angular momentum /, total angular momen- 
tum j = I ± 1/2, and magnetic quantum number m. 

The TH operator rj" of Eqs. ^ and (g]) is obtained 
by keeping the terms with a = a' , 



= X! "5 ' 



(28) 



with Tq, = (a I a) in the TH notation. Therefore, we define 
the correction operator. 



St. 



a~ ba {a'\a) . 



(29) 



Then the correct isospin operator and the exact Fermi 
matrix element are given by 



r 1 = T 



TH 



(30) 

Mp = {f\Tl%) + {f\Sr+\t). (31) 



We evaluate both terms in Eq. (j3ip to second order in 
ISB interactions. This will explicitly demonstrate that 
the second term due to radial excitations is of the same 
(second) order as the TH term. We start with the latter. 



M™ = (/|r™' 



A/o-(/|^at6„(l-r„)|z). (32) 



Here we have applied the TH procedure and neglected 
the isospin-mixing correction 6ci of Eq. (jlOp . which is 
smah in Ref. [1], so that Mj" = Mq for = 1. Since 
(1 — Tq) starts at second order in ISB interactions, we can 
replace the full eigenstates in the second term in Eq. ((32|) 
by \i) and |/). Using a single-particle version of Eq. (PT|) . 
we express (1 — r^) in terms of the matrix elements of the 
one-body potential Uc, that accounts for the difference 
between the |5) and la) basis states, to second order. 



{a'\Uc\a) 



(33) 



where Ea are the eigenvalues of the isospin-symmetric 
single-particle Hamiltonian T+U . For simplicity, we take 
|z) to be a Z — proton plus core configuration. We de- 
fine the occupation probabilities pa of the proton excess 
in the |a) basis, normalized so that pa = (i|r_r_|_|i) = 
2T.^ As a result, we find for the TH term 



M 



TH 



Mn ~ - 



1 1 



2T ^ ' 

a, a.' jta 



\{c^'\Uc\a)\ 

{Ea — Ea' ) 



(36) 



^ To clarify the notation, the symmetry-hmit matrix element can 
then be written as 



Mo 



iwE^s«5E4^,5l*)' (34) 

= 7^wE^5fe3N) = ;^E?- (35) 

where the |5) basis, appropriate for the state was used for 
Tjf. and T_. 



In the limit of sharp occupation probabilities pa, the a 
(a') sum is over occupied (unoccupied) states. 

Next we evaluate the contributions due to radial exci- 
tations, dMp = {f\6T^\i). To second order in ISB inter- 
actions, we have 



SMf « {f\6T+\i) + 
1 



1 



2T 

(i|T_ VcAf 



1 



1 



Ef - AfiJA 



St+ \i) . (37) 



While we will estimate SMp making similar approxima- 
tions as for M™ of Eq. , the result Eq. (|37|) provides 
a correction term that can be included in future numeri- 
cal calculations of ISB corrections. For (5t+ we also need 
the overlap {a'\a) for a' 7^ a, which starts at first order. 



1 



{a'\Uc\a) 



(38) 



Ea — Ea' 

We start with the first term of Eq. ([57)1 . 
{f\6r+\^) = ^(^|^6^a^ ^ al ba {a'\a) \z) , 
i=(z| b\ha{a'\a)\i), (39) 

a^a'^a 



where we have used that, for the considered configuration 
of the neutron annihilation and creation operators 
evaluate to 5~ ^. After transforming ba to the \a) basis, 
using the Hermitian conjugate of Eq. ()25p . we obtain 



2T 
1 



a..Oi'^cx. 



\{a'\Uc\a) 

[En — E„' 



(40) 



We estimate the second and third terms of Eq. pT]) 
using a closure approximation, that is we replace Ei — 
AiHKi by Ai?i < 0, and similarly for Ef - AfHAf. In 
addition, we approximate Vc by the ISB one-body po- 
tential J2-y,j'^-y{l\UcW) b\ bj', where 7 7^ 7' ensures the 
action of the projectors A; and A/. After contracting the 
neutron operators, we find for the second term 



1 



(i| T_ St^- 



1 



■A,Vc\t) 



E, - A,HA 

[a'\Uc\a 



a,a'^a 



|2 



E 



(41) 
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For the third term, we obtain 



1 



2T 
1 



(i|r_ VcA/ 



1 



Ef -AfHA 



dr+ \i) 



2T 



E E {i\UcW)b\b,,b, 



AE 



f 



\^). 



For nonzero overlap, this requires two-particle-two-hole 
configurations in while the first and second terms 
receive contributions at the level of the best Slater de- 
terminant. Assuming residual interactions arc weak, we 
neglect the third term. Combining Eqs. (HO]) and (HI]), 
we have for the correction term 




(42) 



Comparing our estimate 6Mf to the corresponding TH 
term, Eq. (|36p . demonstrates that radial excitations are 
significant. The same estimate, Eq. is found when 

Vc is approximated by the isovector part of the ISB one- 
body potential, J^j^y^-yillUcW) (^ ^' - a^')/2- In 
this case, the second term of SMp is 1/2 of Eq. (|¥T|) . but 
the third term also yields 1/2 of this (for AEi = AEf). 

Assuming radial excitations are dominated by n to 
n-l- 1, we have Ea — Ea' = 'ihio, where w is a typical oscil- 
lator frequency. Moreover, if the relevant excitations are 
dominated by the isovector monopole state, the value of 
\AEi \ ranges between 3 and Ahiu Q. For Ea — Ea' = ihuj 
and \ AEi \ = Afnu, we find that the correction term com- 
pletely cancels the TH contribution Sc2 at second order. 
This result is similar to the energy-degenerate case of 
Auerbach Q. Our estimate shows that, if the contribu- 
tions of the isovector monopole state dominate ISB, the 



radial excitations neglected by TH decrease ISB correc- 
tions. 



V. IMPLICATIONS FOR ISB CORRECTIONS 

We have used the exact formalism of Ref. to calcu- 
late the renormalization of the Fermi matrix element due 
to radial contributions. We expressed the correct isospin 
operator as a sum of the TH operator and a correction 
term involving radial excitations, which were shown to 
be significant and estimated to decrease the ISB correc- 
tions Sc2 of Ref. In addition, for certain conditions 
the ISB corrections due to radial contributions can can- 
cel at second order in ISB interactions. A reduction due 
to the correction term implies that the extracted value 
of Vud nray be reduced. Moreover, our results can pro- 
vide a possible explanation for the smaller ISB correc- 
tions found in Refs. [1, although these calculations 
are more exploratory at this stage. 

An important direction for future research is to include 
the correction term of Eqs. (|29|) and (|37p in numerical cal- 
culations of ISB corrections that follow the TH approach 
or make a similar truncation of basis states. 
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